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p , ■ ABSTRACT 



Or 



OGLE III and MOA-II are discovering 600-1000 Galactic Bulge microlens events 
each year. This stretches the resources available for intensive follow-up monitoring 
of the lightcurves in search of anomalies caused by planets near the lens stars. We 
advocate optimizing microlens planet searches by using an automatic prioritization 
' algorithm based on the planet detection zone area probed by each new data point. 

C/3 , This optimization scheme takes account of the telescope and detector characteristics, 

observing overheads, sky conditions, and the time available for observing on each night. 
The predicted brightness and magnification of each microlens target is estimated by 
fitting to available data points. The optimisation scheme then yields a decision on 
' which targets to observe and which to skip, and a recommended exposure time for 

\& \ each target, designed to maximize the planet detection capability of the observations. 

The optimal strategy maximizes detection of planet anomalies, and must be coupled 
with rapid data reduction to trigger continuous follow-up of anomalies that are thereby 
found. A web interface makes the scheme available for use by human or robotic ob- 
servers at any telescope. We also outline a possible self-organising scheme that may 
be suitable for coordination of microlens observations by a heterogeneous telescope 
0^ . network. 



o 



Key words: gravitational lensing, planetary systems, methods: observational 



o 

> 

X 

H ; 

5r , 1 INTRODUCTION planet anomalies last t p ~ 3(m p /mj) 1/2 days, thus a few 

days for Jupiters or a few hours for Earths. The probabil- 
Gravitational microlcnsing reveals stars and planets that .. , , n „ <w ; a/2 

ity that the planet is detectable is Pd ct ~ 0.2lm„/mi ' 
magnify the light from a background source star (Mao r , , ,„.,,,,, • „ „ r n r> 

„ ° , „ . _ TTT * . for cool planets m the lensing zone , a ~ 0.5 — 2 Re ~ 

& Paczynski 1991). The wide-field OGLE III (Udalski, x _ 4(Mi/Q _ 3 Mq) i/2 atj (GouM & Loeb lg92) 

et al. 2003) and MOA IV surveys of Galactic Bulge 

starfields discover ~ 600 - 1000 microlensing events each When a planet anomaly is well sampled by observations, 

year. During these events, a background source star bright- its duration, timing, and shape determine the mass ratio, 

ens and fades, sometimes by many magnitudes, in t B ~ q = m p /M* and the orbit size a relative to the Einstein ring 

30 (Af*/0.3 M ) 1/2 days as the intervening ~ 0.1-1 M Q radius Re- Roughly speaking, the planet anomaly's duration 

lens star crosses near the line of sight. A planet near the lens t p sets the mass ratio, q ~ (t p /t E f, and its time t, relative 

star acts as a smaller lens, smaller by a factor (m p /M t ) 1/2 . to the event peak at to, measures the projected planet-star 

When appropriately placed, the planet can produce a brief separation a sin /Re ~ (t-t )/t E . The lens star's distance 

but easily detectable flash or dip in the lightcurve. Such D and mass M*, when constrained by the event timescale t E , 

are initially uncertain to factors ~ 3. Several methods using 
finite-source effects, parallax, and proper motion can further 

* http://www.astrouw.edu.pl/~ogle constrain the lens geometry to establish m p , a and M* with 

T http://www.phys.canterbury.ac.nz/moa/ higher accuracy (Gould 2009). For example, high-resolution 
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imaging several years after the event can detect the lens star 
flux, colour and proper motion (Bennett, Andreson, Gaudi 
1996). 

The m 1 / 2 dependence of Einstein ring sizes makes mi- 
crolcnsing more sensitive to low-mass planets than other 
methods. The microlens signature of an Earth-mass planet 
is brief, a few hours, but can be strong enough for easy 
detection (Bennett & Rhie 1996; Dominik, et al. 2007) pro- 
vided one is observing the right star at the right time. With 
a detection probability Pdct ~ 0.2(m p /mj) 1 '' 2 (Gould & 
Loeb 1992), a dedicated survey monitoring ~ 10 3 events 
with < 1 hour sampling could reveal ~ 10 ?7© cool Earths, 
if each lens star has 779 of them. While this level of moni- 
toring has not yet been achieved, significant constraints on 
the abundance of large cool planets were established (Gaudi, 
et al. 2002; Tsapras, et al. 2003; Snodgrass, Home, Tsapras 
2004) even before the first secure microlens planet detection; 
r?j up < 20%. 

A two-stage strategy is currently employed for microlens 
planet searches. The OGLE III and MOA II teams use their 
dedicated wide-angle survey telescopes to discover the mi- 
crolens events. Follow-up teams then deploy networks of 
small narrow-field telescopes distributed in longitude to ob- 
tain more intensive coverage of the most promising of those. 
Two primary strategies are currently advocated and followed 
by the follow-up teams. A strong focus on high magnifica- 
tion events, which have the highest probability of revealing 
planets, is advocated (Griest & Safizadeh 1998; Rattenbury, 

et al. 2002) and put into practice by /jFUN +. The high- 
magnification events are often identified a few days in ad- 
vance, permitting the rapid mobilisation of many telescopes 
to cover the peak of the lightcurve as intensively as possible. 
The PLANET§ Collaboration (Albrow, ct al. 1998) deploys 
a network of small ground-based telescopes to achieve quasi- 
continuous coverage of the most promising events. This ef- 
fort has been joined by RoboNet^ (Burgdorf, et al. 2007; 
Tsapras, et al. 2009), using three 2 m robotic telescopes. A 
much larger robotic telescope network is being laid out by 
LCOGtII in the next few years (Tsapras, et al. 2009). With 
the prospect of this network of 24 0.4m and 18 1.0 m robotic 
telescopes contributing to microlens planet searches, auto- 
mated strategies will be increasingly important to effectively 
organise the follow-up observations. 

The OGLE-2002-BLG-055 lightcurve has one good data 
point that is 0.6 mag high. While this could be a planet 
anomaly (Jaroszynski & Paczynski 2002), undersampling 
prevents adequate characterisation of this event (Gaudi & 
Han 2004). In the first secure characterisation of a microlens 
planet, the lightcurve of OGLE-2003-BLG-235/MOA-2003- 
BLG-053 exhibits two fold caustics separated by 7 days, at- 
tributed to lensing by a ~ 1.6 mj planet (Bond, et al. 2004). 
The 2005 season revealed three microlens planets, ~ 3 mj 
OGLE-2005-BLG-071Lb (Udalski, et al. 2005; Dong, et al. 
2009), ~ 6 m OGLE-2005-BLG-390Lb (Beaulieu, et al. 
2006), and ~ 13 m OGLE-2005-BLG-169Lb (Gould et al. 



t http: //www. astronomy . ohio-state . edu/ ~microf un/ 

§ http://planet.iap.fr 

^ http://robonet.lcogt.net/ 

II Las Cumbres Observatory Global Telescope, http : //lcogt . net 



2006). With two small planets among the first four microlens 
planet discoveries, the abundance of small cool planets must 
be higher than that of the larger cool Jupiters. In 2006, the 
complex lightcurve of the high-magnification event OGLE- 
2006-BLG-109 revealed two planets with mass and orbital 
size ratios that are strikingly similar to Jupiter and Sat- 
urn, scaled to a lower-mass (~ 0.25 M§) host star (Gaudi, 
et al. 2008). From the 2007 season, MOA-2007-BLG-192 
(Bennett, et al. 2008) appears to be a brown dwarf with a 
~ 3 m© planetary companion. Other planets from 2007 are 
not yet published. It appears reasonable on present evidence 
to expect increasing numbers of microlens planet discover- 
ies, leading to detection of cool Earth-mass planets within a 
few years, provided the capabilities for intensive monitoring 
of OGLE III and MOA II events continues to improve. 

This paper develops an optimal strategy for reactive mi- 
crolens planet searches that may help to increase the planet 
discovery rate, particularly with dedicated telescope net- 
works. Section 2 briefly reviews microlens lightcurves to de- 
fine notation and establish a few results for later use. Section 
3 employs numerical integrations and scaling laws to quan- 
tify the detection zone area that we propose as the metric of 
success for a microlens planet search. Section 4 develops the 
optimal observing strategy. Section 5 discusses several prac- 
tical issues, and outlines a possible self-organising scheme 
based on continuously varying target priorities, that may be 
suitable for coordinating microlens observations by a hetero- 
geneous telescope network. Section 6 summarises the main 
results, and describes our web interface to the PLOP (Planet 
Lens Optimisation) algorithm. 



2 MICROLENS LIGHTCURVES 

2.1 Point-Source Point-Lens (PSPL) Lightcurve 

During a microlensing event, light from a background source 
star reaches the Earth along paths that bend toward an 
intervening lens star. With perfect alignment of the observer, 
lens and source, the observer sees the background star as an 
Einstein ring of angular radius 9e centred on the lens. A 
light ray with impact parameter R bends toward the lens 
mass M by a small angle 



25l 
R 



(1) 



where Sl = 2 G Ml /c 2 is the Schwarzschild radius of the 
lens. The point-mass gravitational lens has strong spherical 
aberration, the effective focal length being 



f =R 

J a 



R l 



(2) 



If Dl and Ds are the observer-lens and observer-source dis- 
tances, respectively, the lens formula of geometric optics is 



1 + 1 



1 

7 



2S L 



(3) 



Solving for R gives the radius of the Einstein Ring, 

ife = (2S L L> s X(l-X)) 1/2 , (4) 

where X = Dl/Ds is the lens/source distance ratio, < 
X < 1. The angular radius of the Einstein ring is 



0e 



Re 
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2 5 L Ds 



X 



1/2 



(5) 
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For Ml = M , Dl = 5 kpc, and Ds = 10 kpc, the Einstein 
Ring radius Re ~ 4 AU corresponds to #e ~ 0.8 milli- 
arcseconds. 

With imperfect alignment, one ray on each side of the 
lens reaches the observer. In this case, the lens equation is 
quadratic with two distinct roots, 



Point Mass Lens 



u ± (u + 4) 
u± = ^ '- 



1/2 



(6) 



giving two image positions, u+ = 6+/8e > 1 for the major 
image and u_ = 6-/ 6s, < 1 for the minor image, in terms of 
u = 6/Os, for the unlensed source. Note for future reference 
that 



2 _ T±B 
u ± - — j— 



where T = 

.1/2 



(7) 



+ u 2 _ = u 2 + 2 and B — u 2 , ~ u 2 _ = 



u (w 2 + 4) 

Fig 1 shows the Einstein ring and trajectories of the 
two images during a microlensing event. On the lens plane 
perpendicular to the line of sight, we define cartesian coor- 
dinates x and y with the origin at the lens star, the source 
star moving in the +x direction and crossing the +y axis at 
closest approach. In units of #e, the source- lens separation 
is 



2\l/2 



(2 2 \ J 

u + u x ) 

with Mo the separation at closest approach, and 
., _ l*{t-to) _ (t-tp) 



1i 



(8) 



(9) 



where [i is the relative proper motion, to is the time of closest 
approach, and the event timescale, 



(10) 



is the time to cross the radius of the Einstein ring. 

The image-lens separations satisfy u\u 2 _ = 1, so that, 
as seen in Fig 1, the major image at u + is always outside 
the Einstein ring, while the minor image at u_ remains in- 
side. The major image slides "over the top" of the Einstein 
ring, while the minor image traces a loop inside the ring. 
Both images become brighter as they approach the Einstein 
ring. Each point on the disc of the source star maps to a 
corresponding lensed position on the image. The images are 
thus stretched in azimuth by a factor u± /u and squashed in 
radius by du±/du. With surface brightness conserved, the 
net magnification arising from the increased solid angle is 



A± = 



1 u± 1 


du± 


1 u \ 


du 
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The image magnifications satisfy A + 
magnification is 



(11) 



l + ^4_ , and the total 



.4+ + A. 



u + 



i 2 



T 



u +2 



u (u 2 + 4) 



T72 



(12) 



Since u changes with time, this defines a characteristic point- 
source point- lens (PSPL) lightcurve (Fig 1). Power-law ap- 
proximations (Fig. 2) for large, intermediate, and small u 
are 




x/R E 




1 
x/R E = (t-t )/t E 

Figure 1. Top panel shows snapshots of the unlensed background 
star and of the two distorted images of it that appear on opposite 
sides of the lens during a micro-lensing event. The major image 
passes over the top of the Einstein ring, while the minor image 
executes a loop inside the Einstein ring. The two images become 
compressed in radius but elongated in azimuth as they approach 
the Einstein ring, resulting in a net magnification. This produces 
the characteristic lensing lightcurves shown in the bottom panel, 
where dashed lines show the separate lightcurves of the two im- 
ages, A+(t) and A—(t), and the solid line is their sum, A(t). Note 
that A + (t) = 1 + A_(t). 



2u- 4 u > 2.5 , 
A(u)-1={ u- 2 /3 0.3<u<2.5 
u^ 1 u £ 0.3 . 

Note that A(u) has the inverse: 
^( 2 [(l-^ 2 )- 1/2 -l]) 1/2 . 



(13) 



(14) 



2.2 Finite Source Size 

For a point source at high magnification, as it ^ 0, A ^ 1/u 
becomes formally infinite, corresponding to the formation 
of an Einstein ring of infinite magnification and infitesimal 
thickness. This point-source approximation breaks down, 
however, when the source star's finite size becomes impor- 
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Figure 2. The point-source point-lens (PSPL) magnification 
A(u) compared with the power-law approximations defined in 
Eqn 13. 



tant. In Fig.l, curvature of the highly magnified images is al- 
ready evident. At still higher magnifications, the major and 
minor images extend farther in azimuth, eventually touch- 
ing each other and merging to form an Einstein ring of fi- 
nite width. The magnification remains finite due to the fi- 
nite source size (Bennett & Rhie 1996; Dominik 1998). The 
source star's angular radius 8s = Rs / Ds becomes compara- 
ble to the angular radius 9e of the Einstein ring, at 

= 0s = Es 

(2S L D S X{1 -X)) 1/2 

(15) 

=<"»» (n)terw 5 - 

Finite-source effects set in at high magnification, A £ 1 /us , 
thus (for Ml ~ 0.3 Mq) at A > 500 for a main-sequence 
source star with Rs ~ R®, or already at A > 5 for a giant 
source star with Rs ~ 100 J?©. 

Finite source effects are important for planet anomalies 
when the source star's angular radius 8s exceeds that of the 
planet's Einstein ring (Bennett & Rhie 1996). Since 8 P = 
q 1 / 2 6e, the finite source effect is important when 

-(ft)" (ft)"' (A) - <«) 

Since m© = 3 x 10" 6 M©, this is 

m p = gM L <0.6m©(^) 2 ( T ^ x ) . (17) 

For large source stars the anomaly from a small planet 
can be smeared out and diluted in amplitude, rendering it 
undetectable. On this basis, detection of Earth-mass planets 
is more favourable with main sequence source stars (Bennett 
& Rhie 1996), though a detectable (~ 5%) signal can arise 
even when the source star is a giant (Dominik, et al. 2007), 
provided the planet's alignment with one of the image trajec- 
tories is favourable. We do not consider finite-source effects 
further in this paper. 



2.3 Binary Lens Anomalies 

A planet near the lens star acts like a small defect in the 
gravitational lens. If the planet is well away from the two 
image trajectories, the light it deflects does not reach Earth. 
In this case the planet has no measurable effect on the 
lightcurve and thereby evades detection. However, if the 
planet is close to one of the image trajectories, its grav- 
ity can significantly perturb the bundle of light rays that 
would otherwise reach Earth. This distorts the image and 
changes the magnification to produce a brief anomaly in the 
lightcurve. The magnification curve A2(t) for a star+planet 
lens deviates by a factor 1 + S(t) from the corresponding 
PSPL magnification curve A\(t): 

A 2 (t) = Ai(t)(l + ($(*)) • (18) 

This defines the planet anomaly 8(t), which depends on 
three additional parameters: the mass ratio q, and the co- 
ordinates, x and y, of the planet's projected position on the 
lens plane. 

The planet anomaly may be brief but large. The planet's 
Einstein ring radius is 

r p = R E q 1/2 . (19) 

The planet anomaly may be large when one of the source 
images passes closer to the planet than r p , provided the 
source is not much larger than r p . The duration of the planet 
anomaly is roughly the time it takes the image to cross the 
diameter of the planet's Einstein ring, 

t p =t E q 1/2 . (20) 

Detecting the planet requires data points in the lightcurve 
of sufficient accuracy and at the right time to detect the 
anomaly produced as the image passes by the planet. 



3 PLANET DETECTION ZONES 

3.1 Definition of Detection Zone 

We define the "detection zone" as the region on the lens 
plane (x,y) where the lightcurve anomaly S(t, x, y, q) is large 
enough to be detected or ruled out with high confidence by 
the observations. For N data points with fractional accuracy 
Oi at times U, the detection zone is defined by 

g(te^)) 2 >Ax ', (2D 

for some detection threshold A% 2 ■ This detection threshold 
must be set high enough so that noise affecting the obser- 
vations does not produce false triggers at an unacceptably 
high rate. A\ 2 in the range 25 to 100 corresponds to a 5a 
to 10cr deflection in the lightcurve if the anomaly is confined 
to a single data point. 

Fig. 3 highlights the detection zones for a planet with 
mass ratio q — 10~ 3 derived from a lightcurve A(t) with 
maximum magnification A = 5. Data points uniformly 
spaced in time sample the lightcurve with an accuracy 
a = (5/A 1 ^ 2 )%, and the detection criterion is A\ 2 > 25. 
Each data point probes for planets close to the correspond- 
ing major and minor image positions. If a planet is placed 
inside one of these detection zones, the lightcurve at time 
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Figure 3. Detection zones on the lens plane indicate the regions 
where a planet with mass ratio q = m/M = 10 -3 is detected 
with Ax 2 > 25. The lightcurvc A(t) has maximum magnification 
Ao = 5, and the accuracy of the measurements is a = (5/A 1 / 2 )%. 
Each data point probes for planets close to the two images of 
the background star. The detection zone areas scale roughly as 
U^Rl (2A-I)q/(*A X ). 



U is perturbed by 8 > 5 a = (25/A 1,/2 )%. Improving the 
accuracy of the data or increasing the mass of the planet 
enlarges the size of the detection zone. 



3.2 Numerical Evaluation of Detection Zone 
Areas 

If the data points in the lightcurve are widely spaced, as they 
are in Fig. 3, then the detection zones arising from different 
data points are well isolated from each other. We may then 
evaluate numerically the area O of the detection zone that 
is carved out by each data point. This quantifies the planet 
discovery potential of each data point. 

Fig. 4 shows a close-up of the detection zones defined 
by this criterion. The region displayed is chosen in advance 
from rough estimates and is used for numerical evaluation of 
the dectection zone area Q. The cases shown illustrate how 
the detection zones shrink and change shape as the magnifi- 
cation A(u) declines with increasing lens-source separation 



u. The detection zone shapes are complicated. At small u 
and high A they bear some resemblance to 4-leafed clovers 
with radial and azimuthal lobes straddling the image posi- 
tions. With increasing u, decreasing A, the azimuthal lobes 
of the major image detection zone collapse radially. The ra- 
dial lobes then merge radially to form a circular detection 
zone as u — > oo. On the minor image detection zone, the ra- 
dial lobes merge and vanish, leaving two isolated azimuthal 
lobes that shrink and vanish. 



3.3 Scaling Laws for Detection Zone Areas 

It will be helpful to understand how planet detection zone 
areas scale with the accuracy of the data, the source magni- 
fication, and the mass of the planet. If accurate scaling laws 
can be found, we may then avoid long numerical calculations 
to determine the detection zone area. In this section we de- 
velop a useful analytic formula, and test it against detailed 
numerical integrations. 

Consider first a planet located quite far from the lens 
star, affecting the major image at a time well before or well 
after the stellar lensing event, when u+ » 1 and A+ ~ 1. In 
this case the planet and star act as independent lenses, and 
a significant anomaly occurs when the major image sweeps 
past the position of the planet. If z is the separation between 
the planet and the major image, the planet magnifies the 
major image by a factor A(u p ), where u p — z/r p , and r p = 
q 1 ^ 2 -Re is the planet's Einstein ring radius. 

A data point with fractional uncertainty a can detect 
the anomaly 8 when 

^%)-l>a(A X 2 ) 1/2 • (22) 

This criterion corresponds to a circular detection zone 
around the major image at the time of the observation. The 
radius of the detection zone is found by solving Eqn. (22) 



for 



and hence 



3 u p . Using Eqn. (14) to invert A(u p ), 



the area of the detection zone is 

= (^) 2 =[ 2 ( 1 -( 1 +^r 



7rr: 



-1/2 



- 1 



(23) 



Using the approximations in Eqn. (13), the corresponding 
approximations for the detection zone area are 



irr n 



(2/5) 1/2 5 < 0.05 , 
(3 Sy 1 0.05 < 8 < 3 
S- 2 <5>3. 



(24) 



Fig. 5 indicates that the middle approximation predicts 
fairly accurately the detection zone area for anomalies in the 
range 0.05 < 8 < 3. For Ax 2 = 25 this range corresponds 
to fractional uncertainties 0.01 < a < 0.6, quite appropriate 
for CCD data, giving 

O _ nq 

Mr 37 • 



(25) 



Now move the planet closer to the lens star. The planet 
acts as a defect in the stellar lens, and the resulting detection 
zone can have a quite complicated shape (Fig. 4). We guess 
that the detection zone area may scale roughly with the 
magnification A. However, we find that Q oc 2A — 1 is a 
better approximation: 

Jp^HpA-l). (26) 
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Planet Detection Zones Planet Detection Zones Planet Detection Zones Planet Detection Zones 
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Figure 4. Planet detection zones for mass ration q = 10~ 3 near the major and minor image positions for data points at magnifications 
A = 100, 5, 1.5, and 1.1. The detection zones are symmetric around the image positions in (8,\nu) coordinates. For u — * oo and A — > 1 
the major image detection zone becomes circular and the minor image detection zone splits in two and vanishes. 
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Figure 5. The detection zone area Q scales with the area irr^ 
of the planet's Einstein ring, and decreases with the size of the 
lightcurvc anomaly <5 that can be detected by the data. The ap- 
proximations (dashed lines) are those defined in Eqn. 24. 



Fig. 4 shows that the detection zones are roughly sym- 
metric as a function of \nu and 8 rather than x and y. It 
may therefore be more appropriate to express detection zone 
areas using a dOdlnu metric, evaluating 



J P(detM) dJtJi* , 



rather than 



Re 



P(det|u,0) ududO 



(27) 



P(det\x,y) dxdy . (28) 



The dd d In u metric may be appropriate from a second 
perspective. Exo-planet orbits should have random orienta- 
tions, so the planets should be uniformly distributed in 8. If 



their orbit size distribution is also roughly uniform in log a, 
then the planet distribution on the lens plane will be roughly 
uniform in In u, and the planet detection probability will de- 
cline to zero long after the peak of an event, rather than 
reaching a positive asymptotic value. In effect the dOdlnu 
metric recognizes a detection zone with area Q as more likely 
to include a planet, and therefore more valuable to us, when 
it is measured at small u and probes a larger range of log u. 
For small detection zones the two metrics are related by 



R E u 



7T<7 ( 2A—1 
{—?- 



(29) 



But we must be more careful to treat separately the major 
and minor image detection zones, surrounding the images at 
u+ and u- respectively. This gives 



w ± « || F ± (A) , 

where for the major image at u + 

F+(A) = 7^ - B(T + B) - 
and for the minor image at it_ 



A 



1/2 



1/2 



(30) 



(31) 



(32) 



In deriving the above expressions, we used Eqns. (7), (11) 
and (12) to write 



A± = 



T±B 



T±B 



T = u\ + ul 



2A 



(A* - 1) 



T72 



(33) 



u 2 _ =u (u 2 + 4' 



1/2 



(A 2 - 1) 



T/2 
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Figure 6. The magnification, detection zone area, image velocity, 
and detection zone crossing time are shown for an event with 
peak magnification Aq = 5. The dashed and dash-dot curves are 
for the major and minor images respectively, except for the image 
velocity panel where the dashed curve is dd /dt while the dash-dot 
curve is dlnu/dt for the major image. 




Test of Scaling Law 




Figure 7. The detection zone areas w in the dSdlnu metric 
increase with magnification A, increase with mass ratio q, and 
decrease with the size of the lightcurvc anomaly <5 that can be 
detected by the data. The fast analytic approximation (dotted) 
defined in Eqn. (34) is compared with more exact numerical inte- 
grations. Dotted curves in the lower panel give fractional contri- 
butions of the major and minor images, in the analytic approxi- 
mation. 



The total detection zone area, summing the detection areas 
of both images, is 



7T q 
37 



F(A) 



TYqF(A) 
3 (A X 2 ) 1/2 a(\nA) 



(34) 



with a(lnA) the fractional accuracy in measuring A, A\ 2 
the threshold for planet detection, and 



F(A) =F+(A)+F-(A) 



2T + B 



A-l 



1/2 



(35) 



The functions F±(A), for the separate images and F(A) for 
the total detection zone area are plotted in Fig. 6 for an 
event with Aq = 5. 



3.4 Analytic vs Numerical Detection Zone Areas 

Fig. 7 compares numerically-integrated detection zone areas 
with the analytic result in Eqn. (34) . The figure shows how 
w depends on A— 1 for three mass ratios, q = 1CF 3 , 1CF 4 and 



10 5 , and for three accuracies, S = a(lnA) (Ax 2 ) 1 ^ 2 = 5, 
10 and 20%. 

The analytic approximation clearly captures the main 
scaling, w oc qA/S, for 2 < A < 300. However, it is a 
rough guide rather than a superb approximation. At mod- 
erately high magnifications, 10 < A < 300, the analytic 
and numerical results agree to within ~ 5%. Here we have 
F+ « F- « (2 A— 1)/2, the two images contributing roughly 
equally. At intermediate magnifications, 2 < A < 10, the 
analytic result is low by up to ~ 10%. This A-dependent 
bias could be reduced by adjusting the formula for F(A) in 
Eqn. (35). However, the accuracy is already sufficient for our 
purposes in the range 2 < A < 300. We discuss below the 
breakdown at larger and smller magnifications. 

At very high magnifications, A > 300, our analytic ap- 
proximation breaks down because the detection zones ex- 
tend so far in azimuth that the major and minor image 
zones touch each other and merge together. Since the detec- 
tion zones span a roughly equal range in 9 and in lnu, the 
saturation in 9 should set in when w ~ (2tt q A/38) ;J it 2 , 
i.e. 



3tt 8 



(36) 
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We enter a new regime in which the azimuthally-merged de- 
tection zone may continue to expand in In u but is saturated 
in 9. The slope should drop to w <x (qA/S) 1 '" 2 , i.e. 



Detection Zone Area 



(At I) 



1/2 



(37) 



These expectation are roughly consistent with the behaviour 
in Fig. 7. 

Note that in the very-high magnification regime finite 
source effects will also become important, altering the re- 
lationship between u and A with u > us, as discussed in 
Sec. 2.2. As we neglect both effects in our analysis, our 
scaling law applies only up to a maximum magnification 
A S 8/q ~ 300. 

At low-ish magnifications, A < 2, Fig. 7 shows that the 
analytic formula over-predicts detection zone areas by fac- 
tors of up to ~ 3. The structure is independent of q but de- 
pends on S and A, due to the complicated structure of the de- 
tection zones as seen in Fig. 4. In this regime, A ~ 1 — 2 u~ 4 , 
F+ « -uT 2 , F- « 2u~ 2 , and F w 3u~ 2 . With F- fa 2F+, 
the analytic formula gives the minor image twice the detec- 
tion area of the major image. The F+ formula has correct 
asymptotic behaviour at both high and low magnifications, 
so the problem is with the F- formula. In fact at low mag- 
nifications the radial lobes of the minor image merge and 
disappear, as seen in the right two columns of Fig. 4. This 
cuts the total detection zone area by a factor of about 3 
when A drops below a threshold, 



1 



u 



3<5 . 



(38) 



As shown in Fig. 8, a fairly successful attempt to repair this 
deficit is 



w =M(4ri) (A + (l+A)C(x)) , 
where x w 3 S/(A — 1) « 3 u 4 5/2, and 



C{x) 



0, 



(39) 



(40) 



cuts off the minor image contribution at the appropriate 
threshold. However, this makes F_ depend on 5 as well as 
A. 



3.5 Detection Zone Crossing Timescale 

The cadence of observations aiming to detect planet-like 
anomalies should ideally be matched to the time it takes 
for the images to cross the detection zones, so that each 
measurement probes for planets in an independent region of 
the lens plane, rather than overlapping with the region al- 
ready sampled by previous measurements. This will clearly 
depend on the size of the detection zones, and on the speed 
at which the images move on the lens plane. 

If we approximate the width of the detection zone as 
the square root of its area, then the detection zone crossing 
timescale is 



1/2 



(41) 



where w± — qF±{A)/8 is the detection zone area in the 
A8 d In u metric, and v is the corresponding image speed, 




Test of Scaling Law 




Figure 8. As in Fig. 7 but using the approximation in defined in 
Eqn. (39) rather than Eqn. (34). 



(*)" 



+ 



d In u± 



("■IT 



(42) 



We show below that v is actually the same for both images, 
they move at the same speed on the 9 vs In it plane. The 
image speed and planet detection timescale are plotted in 
Fig. 6 for an event with Ao = 5. 

To evaluate the azimuthal velocity 9, note that both 
images sweep around in 9 at the same rate as the unlensed 
source position. With x = (t — to)/tE and y = un, the un- 
lensed source position moves with 



is 



Differentiate 9 = arctan ^ to find 



da: x + y 



—up 

2 , 

U tE 



(43) 



(44) 



Evaluating the lnw± image velocities is more involved, 
but leads to a simple result. First, differentiate v 2 = x 2 +y 2 
to find 



ax 



x x 
u 



x 
uTe 



Next, differentiate (u±) 2 = (T ± B) /2 to find 



u± 
u± 



T±B _ uu CL±j4) 



4(u±)' 2(u ± Y ' 
where we have differentiated Eqn. (33) to find 
t = 2uii , B = 2uuA . 



(45) 



(46) 



(47) 
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Then, since A = T/B, and (u±) 2 = (T ± B) /2, 



u± 
u± 



uu (B±T) 
B (T±B) 



uu 



x 
UB 



(48) 



Notice that both images have the same velocity in lnu, as 
well as in 8, so that the image velocity v is the same for both 
images. 

Substituting Eqns. (48) and (44) into (42), the image 
speed on the 0, In u plane is 



(vti 



u , x 
u B 



l+4u 2 /u 4 
u + 4 



U +4^0 

u A (u 2 +4) 



(49) 



The image speed is plotted in Fig. 6 for an event with Ao = 
5. The azimuthal velocity 9 dominates near the peak, and 
the radial velocity A\nu/&t dominates in the wings of the 
lightcurve. For large u 2 , the image speed varies as vtE ~ 
1/u. The maximum speed vo £e = 1/ifo is reached at the 
peak of the event, at t — to, where u — uo- 

The planet anomaly crossing timescale is given by 



(fe) - 



w± = q F± {A) u 2 + 4 
{vt E f ~~ 3 1+4^/V 



(50) 



Crossing times for the major and minor images are shown 
in the lower panel of Fig. 6. The crossing time is slightly 
larger for the minor image, by a factor (F- / F+) 1 ' 2 = 
((A + 1) /A) 1/2 (see Eqns. (31) and (32)). The crossing time 
at first rises due to the increasing size of the detection zone, 
and then drops to a minimum at the peak, where the max- 
imum image velocity vo = 1/ (£eUo) is reached. This mini- 
mum crossing time, at the peak of the event, is 

t±(o) =*b«o(§ (a q + 1^1)) 1/2 (^) iIA 

(51) 

-(A)'"- 

Here the final approximation, using F± « A and u ?s A -1 , 
holds to 15% or better for Ao ~ 3. 

For a specific example, consider the crossing time for 
a Jupiter- like planet with q — 10~ 3 , in a typical event 
with £e = 30 d. For a peak magnification Ao — 5, as 
in Fig. 6, we have u = 0.2, F+ = 5 (4/6) 1/2 = 4.1 
and F- = 6 (4/6) 1/2 = 4.9. For a good data point with 
a (In A) = 1%, and a detection threshold at Ax = 100, 
the smallest detectable planet-like anomaly deviates by S = 

a(ha.A) (Ax 2 ) 1//2 = 0.1. The crossing time for the major 



image is 



(52) 



t+(0) = 30 x 0.2 (^$^5) (|) =1.2d. 

For an Earth-mass planet and a typical lens mass Ml ~ 
0.3 Mq the mass ratio is q — 10~ 5 , and the crossing time is 

M0) - 2.3 h (j&j) (^) V2 (£) ~ 1/2 (4) ~ 1/2 .(53) 



4 OPTIMIZING A MICROLENS PLANET 
SEARCH 

4.1 The Observer's Dilemma 

In this section we consider how an observer might try to 
optimize a microlens planet search. We assume that the 



observer has many targets to choose from. This is a good 
assumption because MOA II and OGLE III are finding 
~ 600 — 1000 events each year. The observing time available 
on each night during the winter months when the Galactic 
Bulge is visible from a southern hemisphere site is of order 
10 hours. A dedicated agile telescope spending 2 minutes 
per target could in principle visit over 100 targets per night. 
However, why should equal time be devoted to all targets? 
Surely the brighter and higher-magnification targets warrant 
more attention. By skipping fainter and/or weakly magni- 
fied sources, we can spend more time on the more favourable 
ones. The resulting planet detection zones, accumulated over 
all targets during the night, will then be larger, increasing 
the chances of discovering a planet. At the other extreme, 
when one source is very highly magnified, should we attend 
exclusively to that source, and skip all the others, or should 
we reserve some time for a few of the other sources as well? 
This is the microlens observer's perpetual dilemma. The so- 
lution we propose is to observe always in a way that aims 
to maximise the probability of planet discovery. 

4.2 Accuracy of Photometry 

Because detection zone areas scale as w oc a -1 , e.g. 
Eqn. (34), a critical issue is the accuracy of photometric 
measurements that can be achieved, and the rate at which 
that accuracy improves with exposure time. We assume that 
the data analysis is close to optimal, so that photon count- 
ing statistics dominate the noise budget. Thus CCD readout 
noise, cosmic ray hits, and other noise sources are neglected 
in comparison with the Poisson noise from detected star and 
sky photons. The signal-to-noise ratio then increases as the 
square-root of the exposure time, 

= (si) V2 > ( 54 ) 

where At is the exposure time, and r is the exposure time 
required to reach a signal-to-noise ratio of 1. 

The parameter r controls the exposure time needed to 
obtain information on the current magnification A. It de- 
pends on the telescope collecting area, the detector sensi- 
tivity and bandwidth, on the brightness of the magnified 
source star, and the degree of dilution of its photons by sky 
background and by other stars that are blended with it. In- 
cluding these three sources of Poisson noise, 

f* + /B + /sky 



(55) 



where /*, /b and / s k y are the number of detected photons 
per unit time from the magnified source star, from the lens 
star and other stars blended with the source star, and from 
the sky, respectively. We elaborate these three Poisson noise 
sources below. 

For a star of magnitude m* and spectral energy distri- 
bution f\ (A), we observe thru the atmospere with transmis- 
sion T(A), with a detector effective area A e g(X). The photon 
detection rate can be evaluated precisely as 



f A (A) dA 
J h » 



A eff (A) T(A) , 



or approximately as 
/*«/o(A)T(A)10- a4m * 



(56) 



(57) 
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The photon detection rate from Vega (magnitude 0) is 
/„(/)« 500s" 1 (^) , (58) 

for a telescope with mean effective area A cS over a band- 
width AA near the I band, where most microlens observa- 
tions are taken (for the V band, Vega's flux is 1000 rather 
than 500 photons cm _2 A _1 s _1 ). 

For the source star, magnified by a factor A, 



/* = f s at = forlo- 



rn 



Stars blended with the magnified source star contribute 
Poisson noise to the measurement. The source flux fs and 
blend flux /b are normally measured by fitting observed 
lightcurves (corrected for atmospheric transmission) with 
the model 



f(t) = fs A{t) + /b 



(60) 



When using differential flux measurements Af(t), obtained 
by a difference image analysis, the reference flux added to 
these is somewhat arbitrary. As a consequence, the blend 
flux /b arising from the lightcurve fit is also somewhat ar- 
bitrary, and can even be negative in some cases. 

The blend flux contributing to the Poisson noise in- 
cludes not only flux from the lens star, and any other stars 
that are "exactly" coincident on the sky with the magnified 
source star, but also stars that are close enough on the sky 
so that the point-spread functions overlap. If rrii is the mag- 
nitude and 6i is the angular separation of star i from the 
target star, the blend flux contributing Poisson noise to the 
measurement is 



/ B (A)=/,]Tl0- 



4 (mj-m 



.) e -(6i/Ay 



(61) 



This expression assumes a gaussian point-spread function 
with standard deviation A, and optimal extraction to mea- 
sure the target star flux. Note that the Poisson noise due to 
blended stars increases with the seeing. Although it is not 
yet done in practice, the specific dependence on seeing for 
each microlens target can be evaluated in advance from a 
good-seeing image of the starfield, for example the OGLE 
or MOA finding-chart images made available for each event. 

Finally, the detection rate of sky background photons 
overlapping with the target star is 



/ sky = f AO 2 10-°' 4 



(62) 



where ^i s k y is the magnitude of a square arcsecond of sky, 
and A9 2 is the solid angle subtended by the photometric 
aperture (for aperture photometry) or by the point-spread 
function of the star images (for psf-fitting photometry) in 
square arcseconds. The sky brightness, including e.g. air- 
glow, zodiacal light and scattered moonlight, may be eval- 
uated using a sky model, e.g. (Krisciunas & Schaefer 1991; 
Patat 2003). The effective sky coverage of a gaussian point- 
spread function is 



(63) 



A * 2 = 4 - A2 =2Tn2 W '> 
where A is the standard deviation and W is the full- width at 
half-maximum (FWHM) of the gaussian point-spread func- 
tion. Atmospheric seeing is usually reported in terms of W. 

Combining the above equations, we can rewrite 
Eqn. (55) as 



T* + TB + T s ky , 



with 



10' 



0.4m, -9. 7 



4(»i-">.) c -(9,/A)- ! 



and 



4nA 2 



10 



-0.4( 



Msky" 



(64) 
(65) 
(66) 

(67) 



r sky - T+ -jrjjy 

These expressions make explicit how r depends on the mag- 
nified source star brightness (magnitude m*), on the nearby 
stars blended with the target (magnitude rrn, separation 
6i), on capabilities of the telescope (effective area A cS , 
bandwidth AA), and on observing conditions (sky bright- 
ness Psky, seeing A, atmospheric transmission T). When the 
sky and blend fluxes are negligible, a 100s exposure with 
A e g = 1 m 2 and AA = 10 3 A reaches 1% accuracy at mag- 
nitude m, = 16.8. 

4.3 Optimal Exposure Times 

The "worth" of an observation, from the perspective of 
planet hunting, is proportional to the area of the result- 
ing detection zone. Combining Eqns. (34) and (54), we see 
that detection zone areas increase with the square-root of 
the exposure time, 



At t 
tAx 



1/2 



g {At) 1 ' 2 . (68) 



The proportionality constant g characterizes the "goodness" 
of observing this particular target, 

,. _ q F(A) 



(r A X 



2\l/2 



(69) 



These g values can be used to prioritise the events that are 
available at any given time. They depend on the properties 
of the event, characteristics of the telescope, and on the 
present observing conditions. 

The key point to note here is that the fractional mea- 
surement error decreases with the exposure time, a oc 
At -1 / 2 , and this expands the detection zone area as w oc 
At 1 ' 2 . The detection zone area grows most rapidly at the 
beginning of the exposure, with diminishing returns as the 
exposure progresses. For this reason at some point it be- 
comes advantageous to abandon observations of this target 
in favour of moving on to a fresh target that has not yet 
been observed. 

Suppose that we are contemplating making observations 
of N targets during an upcoming night in which we expect to 
have available a total observing time t. How much exposure 
should we devote to each target? For each target i we can 
calculate the goodness factor gi. If we observe target i with 
exposure time Ati, then 



JV 

t = ^ Au 



(70) 



is the total observing time. The total worth of observing the 
N targets is 
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Dividing time between 2 targets 




At 



t 



At. 



Figure 9. When observing N microlens events, for which the 
planet detection zone areas tu, grow with exposure time At; as 
Wi = gi At; 1 / 2 , the total detection zone area W = . Wi is max- 
imised when the available exposure time is divided in proportion 
to g?. Optimisation for the case N = 2 is illustrated here. 



N 



Wn = (At,) 1/2 



(71) 



Given a fixed total observing time t, we can optimize 
the exposure times by solving 8Wn /dAti = 0. For example, 
with N = 2 targets, as illustrated in Fig. 9, the total time 
is t = Aii + At2 , and the sum of the detection zone areas is 



Wn = gi (Aii) 1/2 + g 2 (t - Ati) 1/2 
Maximizing Wn gives 

8W N _ gi _ ff2 



= 



2 (Aii) 1/2 2 (At 2 ) 1/2 



(72) 



(73) 



and thus Ati/At2 = (31/^2) . Similarly, for the general case 
of iV targets, the optimal exposure times that maximize Wn 
also satisfy Ati oc gl, and are therefore given by 



At 
where 



(74) 



(75) 



If we adopt the optimal exposure times, substituting 
Eqn. (74) into (71) gives the total worth of the observations 



W N = G N t 1/2 



(76) 



This analysis suggests that the optimal strategy to max- 
imize the planet detection capability is to observe all avail- 
able targets, spending more time on the best targets, using 
exposure times proportional to the square of the goodness, 
Ati oc gf. The optimal observer skips no targets. As we will 
see, however, this conclusion is altered when we take account 
of observing overheads. 





time (s) 




time (s) 



time (s) 



Figure 10. Illustration of the reduction in the net exposure time 
and the corresponding degradation in planet hunting capabil- 
ity for a 100s observation accounting for a CCD readout time 
^rcad = 10s, for splitting the exposure into 3 sub-exposures to 
avoid saturation, and for a telescope slew time t s i ew = 30s. 



4.4 Effect of Overheads 

In practice the CCD camera takes a finite time tread to read 
out, and the telescope takes a finite time t s i ew to slew from 
one target and settle into position on the next. Typical read- 
out and slew times are tread ~ 10 — 20 s and t s i ew ~ 1 — 3 min. 
To avoid CCD saturation, a single long exposure may need 
to be broken up into a series of n shorter exposures. These 
overheads reduce the on-target exposure time accumulated 
during an observation time t to 



At = t — t s]ew — n t, 



cad 



(77) 



These overheads diminish the planet hunting capability of 
the observations, as illustrated in Fig. 10. We must allow for 
these overheads when implementing an optimal observing 
strategy. 

If the CCD exposure is too long, the target will saturate. 
If the CCD exposure is too short, the readout noise will 
dominate over sky noise and information will be lost. These 
considerations set the range that should be considered for 
the CCD exposure time: 



train ^ ^cxp ^ til 



(78) 



A total exposure longer than £ max is accumulated by taking 
a series of n shorter exposures, where 



^max H~ ^rcad 



< n < 



t ^slcw 
^min ~\~ ircad 



(79) 



For bright targets where t max , the longest exposure that 
avoids saturation, is less than t m in, the shortest exposure 
that avoids readout noise domination, the need to avoid sat- 
uration must take precedence. Having n > 1 protects against 
cosmic ray hits. Once n is decided, the duration of each ex- 
posure is 



^cxp — it t s \ cw ) /n tread 



(80) 



Allowing for these overheads, the detection zone area be- 
comes 



^ ( ^cxp ~F 



1/2 

2iP (t 

tread ' V 



islew) 



1/2 



(81) 
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Dividing Time among N Events 




20 40 60 

N = Number of targets 



100 



Figure 11. The number of microlens events to observe is taken 
to maximise the total worth of the observations. With too many 
targets, observing time is reduced by the time required to slew 
the telescope from target to target. This optimisation of N is il- 
lustrated for 100 targets with an exponential distribution of good- 
nesses, and for different total amounts of telescope time. With less 
time fewer events are observed, and the total worth is reduced. 



The first bracket accounts for the reduction in on-target ob- 
serving time due to the CCD readout time. We can absorb 
this term into the definition of g, 

n _^ 9 



H72 



(82) 



(1 + (t rca d Acxp)) 

as shown by the dashed curves in Fig. 10. The effect is to 
suppress interest in observing targets that are so bright as 
to require inefficient observations with t cxp < t roa d- A target 
too bright for efficient observations with a large telescope 
may thus remain a prime target for smaller telescopes. In 
this way the scheme may serve well to coordinate observa- 
tions by a community with a variety of telescope types. 

The second bracket in Eqn. (81), allowing for the slew 
time, delays the onset of detection zone growth while the 
telescope is moving from one target to the next. We will see 
below that this term dictates which of the less promising 
targets to omit from the observing schedule. 

4.5 Dividing Time among N Events 

If we try to observe too many of the ongoing events, we will 
spend all night slewing from target to target and no time 
at all collecting photons from the targets. If the total time 
available for observations is t, and slew time is t s i ew , then the 
maximum number of targets we can contemplate observing 
is 



N n 



t/tslew 



(83) 

targets, the total worth of the 



If we observe N < iVmax 
observations will be 

W N = G N (t-Nt slcv ) 1/2 . (84) 

We would like to maximize Wn- The first term Gjv increases 
with N, and the second term (t — Ntsicw) 1 ^ 2 decreases with 
N. Therefore Wn has a maximum value for some N < iV max . 
This is the number of targets that we should observe to 
maximize the planet hunting capability of our observations. 



To make Wn grow as fast as possible, sort the targets 
and consider them in order of decreasing goodness, qn > 
9n+i- We should keep target N + 1 only if Wn+i > Wn- To 
decide whether or not to retain target N + 1, note that 



= 1 + 



(Wn+i 
\~Wn~ 



Target N + 1 survives only if 



^slcw 
iV t s l cw / 



9N+1 > 



Gn t 



1/2 



(t-Ntslw) 



1/2 



(Nn 



N) 



1/2 



(85) 



(86) 



where we have used iV max = t/t B \ evr . 

To illustrate this optimisation, Fig. 11 shows the result 
of dividing time among N = 100 targets with an exponential 
distribution of goodnesses gi, for several different total avail- 
able observing times t, corresponding to iV max = 20,40, .... 
As available time t increases, the optimal strategy spends 
more time on each target, and also extends time to addi- 
tional lower-priority targets. 

To further illustrate, more realistically, we consider in 
Fig. 12 the recommended observations from among 443 
OGLE events that were available on 2003 Aug 31. To de- 
cide on the observing strategy, we first fit a PSPL lightcurve 
model to the OGLE data on each event to evaluate the event 
parameters. This results in predicted magnitudes and mag- 
nifications for each target on each night in question. We 
next evaluate the goodness factors gi for a telescope with 
effective area A e g — 1.1m 2 , with a sky magnitude 19. We 
assume a slew time i s i ow = 60s, a readout time i rca d = 10s, 
a maximum exposure time t cxp < 600s. 

In the top panels of Fig. 12, the dashed curve shows 
how Wn would increase monotonically if there were no slew 
time. The solid curve shows show how Wn at first increases 
with N but then decreases as slew time becomes important. 
For t — 1.5h (left panel of Fig. 12), the maximum number 
of targets that could be observed is iV max = t/t B \ ew = 90, 
but the optimal sampling to maximise Wn undertakes ob- 
servations of just N = 6. For t = 8.5h (right panel) all 443 
targets can be observed, but the maximum of Wn occurs at 
N = 20. 

In the bottom panels of Fig. 12, the resulting lightcurves 
are shown when this strategy is employed on every night. 
The area of the plot symbols are proportional to the observ- 
ing time allocated to each target. On most nights the opti- 
mal sampling spreads observing time over many targets. On 
a few nights when one very high magnification event is avail- 
able, that target captures most or all of the recommended 
observing time. One bright target receives some attention 
even though its magnification is small. The observations in- 
clude fainter targets when more observing time is available. 
Targets fainter than the sky are seldom scheduled. 



5 DISCUSSION 

5.1 Detection is not Characterisation 

We must emphasize that the scheme outlined above is de- 
signed to detect anomalies, not to characterise them. Ob- 
serving many targets for the recommended exposure time 
can be advocated only so long as each new observation in- 
dicates that no significant anomaly is underway. It is there- 
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Figure 12. Optimal sampling of 443 OGLE events available on 2003 Aug 31. Exposure times are chosen to maximize the total planet 
detection zone area, for a 2 m telescope with fixed total observing time 1.5h (left) and 8.5h (right) per night. Observing too many targets 
is inefficient because of the 120s telescope slew time. Observing too few targets is inefficient because planet detection zone areas grow 
only as t 1 / 2 . On the resulting lightcurves (lower panels), the plot symbol areas are proportional to the allocated exposure time. On most 
nights observing time spreads over many targets. On some nights one high-magnificaiton target captures most or all of the attention. 



fore best if a rapid reduction of each new observation can be 
undertaken with sufficient accuracy and reliability to check 
each new data point for consistency or otherwise with the 
PSPL model. This is feasible because only one or at most a 
few stars on each CCD image will be undergoing microlens- 
ing at a given time. The sub-image around the target of in- 
terest can be quickly reduced to measure its brightness with 
respect to nearby comparison stars. In practice the real-time 
image-subtraction pipelines currently in use by PLANET 
and RoboNet can reduce each CCD image within a few min- 
utes of the end of the exposure. 



An implementation of this is the SIGNALMEN anomaly de- 
tector (Dominik, et al. 2007). 

By following this two-stage approach - prioritised multi- 
target anomaly hunting punctuated by episodes of anomaly 
confirmation and characterisation - we can simultaneously 
maximize the opportunity to detect anomalies by observing 
a large number of targets, while retaining the ability to re- 
liably establish the nature of anomalies that we detect. If 
the second-stage is omitted from the observing strategy, the 
risk is a series of single-point anomalies will be found whose 
identity cannot be securely established. 



Whenever a significant anomaly is identified, the ob- 
server can temporarily suspend the anomaly-hunting strat- 
egy of observing many targets in sequence, returning to the 
target that offered up the anomalous data point. Additional 
observations of this target then aim to establish either that 
an anomaly is in progress, or else to dismiss the false alarm 
caused by unreliable data. If the return observations fail 
to confirm the anomaly, then the anomaly-hunting obser- 
vations can resume. If the return observations confirm the 
anomaly, then continuous observations are initiated to clar- 
ify the nature of the anomaly, and an alert can be issued to 
trigger follow-up observations on other available telescopes. 



5.2 Dynamic Priorities 

When event parameters change significantly during a night, 
or when the image positions fall inside detections zones from 
previous observations, then the results derived in the previ- 
ous section are no longer strictly valid. Movement of the im- 
ages will increase, while overlap will decrease the detection 
zone areas. We are nevertheless hopeful that the optimisa- 
tion scheme advocated above will still be helpful in guiding 
follow-up observations. 

One way to cope with the more general situation is 
to employ a scheme with continually-evolving target prior- 
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Figure 13. Top panel shows 5 data points at times ti with accu- 
racy <7i = 5% on the decline of a lensing lightcurve with maximum 
magnification Aq = 20. Bottom panel shows the evolving prior- 
ity given to proposed new data points with accuracy a = 1, 2, 
and 4% (top to bottom) when searching for planets with mass 
ratio q = 10~ 3 (solid curves) and 10~ 4 (dashed curves). The 
priority for an isolated data point is proportional to the planet 
detection zone area, Q oc q(2A(t) — l)/cr, where A(t) the time- 
dependent magnification. The priority drops when the new ob- 
servation would probe for planets inside a detection zone already 
established by previous observations. The priority then recovers 
on the timescalc needed for images to cross the detection zones. 



ities. The highest-priority target is observed. The priority 
of that target must then fall dramatically, since immediate 
re-observation would probe for planets inside the detection 
zone just carved out. The priority should then recover in due 
course, as changes in the event geometry move the image po- 
sition outside of the detection zone. Such a scheme may be 
ideal for fully automatic follow-up observations with robotic 
telescopes, but could also be used by human observers will- 
ing to follow directions from a computer programme. 

A dynamical priority scheme of this sort is illustrated 
in Fig. 13. The top panel of Fig. 13 shows a set of data 
points with accuracy Oi = 5% at times ti during the decline 
of an event with peak magnification Aq — 20. The lower 
panel shows the evolving priority given to a proposed new 
data point at time t with accuracy a = 1, 2, and 4% when 
searching for planets with mass ratio q = 1CF 3 and 1CT 4 . 



The priority is evaluated numerically as the increase in de- 
tection zone area arising from the proposed new data point. 

The dips in priority evident in the lower panel of Fig. 13 
indicate the reduced planet hunting capability caused by the 
overlap of detection zones when the new data point probes 
for planets inside the detection zone of an earlier measure- 
ment. We see in Fig. 13 that the reduction is small for 
a — 1% and substantial for a = 4%. This is because the old 
a = 5% data are important when the new data point is of 
similar accuracy, but unimportant when the new data point 
has much higher accuracy. We see also in Fig. 13 that the 
priority recovery time is faster for q — 1CP 4 than for 1CP 3 . 
This is plausible since detection zone sizes scale as q 1 ■ The 
rather irregular recovery arises from the complicated shapes 
of the detection zones (Fig. 4) . 

In practise there will be not just a single previous mea- 
surement with accuracy a, but rather a set of prior measure- 
ments at times ti with accuracies Oi. Noting that indepen- 
dent measurements combine optimally with 1/a weights, 
the net effect at time t of all prior measurements may be 
approximated by using the scheme 



1 



E 



M [(t - U) [sj\ 



(87) 



where Si is an "expiration time" for the observation at time 
ti, and M(x) is a "memory function", 1 for t — ti and de- 
creasing to for t >> ti, effectively forgetting sufficiently 
old observations. Possibilities for the memory function are 
Gaussian or Lorentzian: 



M(x) 



-x 2 /2 



M(x) 



1 + x 

The detection zone crossing, time worked out in Section. 3.5, 
provides a suitable expiration time s;. 

The new detection zone area grows more slowly due 
to overlap with earlier zones. It is as if an exposure time 
t donc has already been done to achieve the accuracy o(t) = 
(T/idonc) 1 ^ 2 - The new exposure time At then adds to tdonc, 
increasing the detection zone area by 
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Another relevant consideration is that slew times are 
not equal for all targets. The slew time is zero for the cur- 
rent target, and for other targets should increase with their 
angular distance from the current target. If we include the 
slew time, then the increase in detection zone area is 

Aw = g((ma. X [0,At-t slcv ] + t donc ) 1/2 -tli 2 ^ . (91) 

If we require the exposure time to be not shorter than some 
minimum time, perhaps some multiple of the CCD readout 
time, in order to have a reasonably high observing efficiency, 
then one compares the options of observing longer on the 
present target without slew time vs slewing to another tar- 
get. As tdone increases on the current target, the potential for 
increasing its detection zone area declines until it becomes 
better to slew to and expose on the next target. 

This scenario is illustrated in Fig. 14. Here three targets 
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Figure 14. The slew time is zero for the current target, but sig- 
nificant for two others. Observations of the current target should 
continue until the marginal improvement in detection zone area 
becomes less for this target than for one of the alternatives, taking 
the slew times into account. 



are considered. We are currently exposing on target 1, with 
slew times of 100s and 200s to reach targets 2 and 3. We 
consider a minimum exposure time of 40s. The circles show 
the result of slewing to an alternative target and exposing for 
the minimum time. In the top panel we have accumulated a 
1000s exposure on target 1. The circles for both alternative 
targets are below the solid curve, so we should not slew. In 
the bottom panel we have accumulated a 2000s exposure on 
target 1, and this increase in tdone reduces the slope of the 
solid curve to such an extent that the circle on target 3 is 
now just above it. At this point we should therefore decide 
to slew and expose on target 3, rather than remaining on 
target 1. This cycle may be iterated throughout the night. 

We expect an optimisation scheme based on approxi- 
mations like those described above to be helpful in deciding 
how long to continue observing the present target, and which 
target is the best one to observe next. We have not yet simu- 
lated this possibility in great detail, but outline the concept 
here as a possible starting point for a self-organising scheme 
that may be suitable for coordinating optimal microlens ob- 
servations by a heterogeneous network of telescopes. Assum- 
ing rapid sharing of information among the telescope nodes, 
each telescope can independently decide which target to ob- 



serve next, taking into account prior observations made by 
all other telescopes, with their various times and accuracies. 



5.3 Uncertain Event Parameters 

The event parameters to, £e, and Ao are often uncertain and 
correlated in the early stages of an event before the observa- 
tions have sampled both sides of the lightcurve peak. With 
highly uncertain event parameters, large errors may arise 
in the assigned target priorities. A frequent example occurs 
when an early fit to the rising part of the lightcurve suggests 
a very high magnification event that later turns out to be 
of only modest magnification. How will such uncertainties 
affect our strategy? 

One happy aspect: the detection zone areas depend on 
current values of the magnification A and star magnitude 
m*, rather than on the values at the peak of the lensing 
event. This is helpful in the early stages of an event when 
the eventual peak magnification is still difficult to predict. 
On the other hand, in a highly-blended event the true mag- 
nification of the source star can be higher than the apparent 
magnification. 

The event parameter uncertainties remaining after fit- 
ting the PSPL lightcurve model to the extant data points 
can be quantified, for example by using the parameter co- 
variance matrix or Markov-Chain Monte-Carlo techniques. 
The corresponding uncertainy in the event priority may then 
be taken into account using a Bayesian average over the pos- 
terior probability distributions. 

One may also contemplate giving priority to observa- 
tions aiming to reduce uncertainty in the event parame- 
ters. This secondary goal will then need to be traded-off 
in some satisfactory way with the primary goal of discov- 
ering planets. An anomaly found on the rise should attract 
attention, making it likely that accurate event parameters 
will be nailed down by observations across and after the 
peak. For an event well past the peak, it may be too late for 
additional observations to pin down uncertain event param- 
eters, or additional observations at critical stages may help a 
lot to break the ambiguity between blending, magnification, 
and event timescale. Targets could be given reduced prior- 
ity when their event parameters are uncertain and there is 
little prospect of improving them, or higher priority when a 
critical observation would help to nail down the uncertain 
prameters. This issue needs careful investigation. 

Our estimates of detection zone areas assume that the 
underlying lens parameters are or will be well constrained by 
observations outside the planet anomaly. When this is not 
the case, then the actual detection zones will be smaller, be- 
cause the loose event parameters can shift the model toward 
the anomalous data points that would otherwise be able 
to detect or rule out planets. Fig. 15 illustrates this effect, 
where the reduction of detection zone areas is considered 
for the OGLE data on OGLE-2005-BLG-390. In this event, 
one OGLE data point occurs during a planet anomaly. The 
reduction of detection zone areas is noticeable but not large 
enough in cases such as this to be a serious problem for our 
optimisation scheme. It would be a more serious problem for 
events with only a few measurements covering the magnified 
part of the lightcurve. 
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test390 A„=3.1 ± 0.027 P„=4.6 P = 0.00016 



5.4 Targeting Specific Types of Stars and Planets 
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Figure 15. Top panel: OGLE III observations of OGLE-2005- 
BLG-390. One data point occurs during a planet anomaly. The 
point-source point-lens (PSPL) model fits 5 parameters, the peak 
time to , peak magnification Aq , Einstein radius crossing time *e , 
the source flux /g and blend flux /b . Middle panel: Greyscale rep- 
resentation of X 2 ( x > y)> moving a planet with mass ratio q = lO -3 
on the (x, y) lens plane, holding fixed the 5 PSPL parameters. The 
X 2 increases by 100 or more in the white areas, where the planet is 
ruled out, and decreases by 100 or more in the black areas. Bot- 
tom panel: The x 2 (x,y) map re-fitting the 5 PSPL parameters 
for each planet position, showing the smaller size of the resulting 
planet detection zones. 



As our knowledge of the exo-planet distribution function 
accumulates, one might contemplate introducing a prior on 
the parameters q, M*, and a in order to target the planet 
search toward particular types of stars or planets. For exam- 
ple, since ts oc M^ 2 , fast events correspond on average to 
lower-mass stars. Similarly, since a oc tsu± , larger orbits can 
be targeted by observing slower events and observing longer 
after the event peak. It is straightforward to tilt the search 
toward any specific parts of parameter space. However, at 
this stage our knowledge of the cool planet distribution is 
so scant that it is probably premature to invest much effort 
into such fine-tuning. 



6 SUMMARY 

OGLE III and MOA II are discovering 600-1000 Galactic 
Bulge microlens events each year. This stretches the re- 
sources available for intensive follow-up monitoring of the 
lightcurves in search of planets near the lens stars. We ad- 
vocate optimizing microlens planet searches by using an au- 
tomatic prioritization algorithm based on the planet detec- 
tion zone area probed by each new data point. We evaluate 
detection zone areas numerically and validate a plausible 
scaling law useful for rough but rapid calculations. The pro- 
posed optimization scheme takes account of the telescope 
and detector characteristics, CCD saturation, readout time, 
and telescope slew time, sky brightness and seeing, past ob- 
servations of microlensing events underway, and the time 
available for observing on each night. The current brightness 
and magnification of each target are estimated by extrapo- 
lating fits to previous data points. The optimal observing 
strategy then provides a recommendation of which targets 
to observe and which to skip, and a recommended exposure 
time for each target, designed to maximize the planet de- 
tection capability of the observations. This must be coupled 
with rapid data reduction to trigger continuous follow-up 
observations whenever an anomaly is detected. It is hoped 
that the algorithm will provide helpful guidance to follow- 
up observing teams, and may be a useful starting point for 
optimising fully-robotic microlens planet searches. 



6.1 WEB-PLOP 

An implementation of this optimisation scheme, Planet 
Lens Optimisation (PLOP or web-PLOP), can be found at 
http://www.artemis-uk.org/web-PLQP/ (Snodgrass, et al. 
2008). This system was designed with two motivations: to 
provide an optimal target list for the automated observing 
of the RoboNet project (Burgdorf, et al. 2007; Tsapras, et 
al. 2009), and also to provide such lists to human observers 
at any telescope. It is formed of two parts. First, a user 
interface web form takes input for the telescope and observ- 
ing conditions parameters (A c s, isicw, ^rcad, Msky, A(? etc.) 
and the total available observing time, t. Secondly, a back- 
ground code keeps track of the current data on each event 
(from OGLE, MOA, RoboNet, and all teams that make data 
available in real time), and produces a new PSPL fit when- 
ever new data arrives. The results from these fits give the 
event parameters Ie, Aq etc. that are used to predict the 



© 0000 RAS, MNRAS 000, 000-000 



Optimising Microlens Planet Searches 17 



magnification at the requested time of observation. These This paper has been produced using the Royal Astronomical 
two sets of inputs allow the calculation of gt for each event Society/Blackwell Science P^TjrjX style file, 
using Eqn. (69), and therefore an optimal list of targets with 
suggested exposure times for the requested telescope at the 
requested time. The list is then put out in either a machine 
readable or sortable human friendly format. With RoboNet, 
this output controls the telescope, and new data is fed back 
into the PSPL model to close the loop and give priorities 
that are based on data just taken. For human observers at 
other sites, the output pages are customisable to display 
any desired parameters along with the priority of each mi- 
crolensing event, and also show light-curves and detection 
zone maps along with links to the finding charts and orig- 
inal OGLE and/or MOA pages for each. Although written 
for RoboNet, this prioritisation tool is freely available and 
other microlensing observers are encouraged to make use of 
it. 
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